arXiv:1503.04872v2 [astro-ph.CO] 7 Oct 2015 


Transients in finite infiation 


Andrew Scacco and Andreas Albrecht 
Department of Physics, UC Davis 
(Dated: October 9, 2015) 

We test a model of inflation with a fast-rolling kinetic-dominated initial condition against data 
from Planck using Markov chain Monte Carlo parameter estimation. We test both an potential 

and the R + E? gravity model and perform a full numerical calculation of both the scalar and tensor 
primordial power spectra. We And a slight (though not significant) improvement in fit for this model 
over the standard eternal slow-roll case. 


I. INTRODUCTION 


One of the greatest sources of data for modern cos¬ 
mology is the cosmic microwave background radiation 
(CMB). This has been measured to extreme precision 
and the concordance model of cosmology has achieved 
tremendous success in matching the data. However, there 
are still small anomalies, one of which is a deficit of power 
in the CMB at low multipoles. A cutoff in the primor¬ 
dial power spectrum from inflation at small k translates 
into a reduction of power in the CMB at low multipoles. 
This low-.^ anomaly consists mainly of a slight dip in the 
power spectrum at € ~ 20 — 40 and a slightly reduced 
quadrupole, though because of cosmic variance the sig¬ 
nificance of the anomaly is not high. In light of this low-£ 
anomaly, there has been a great deal of interest in the¬ 
ories that predict a cutoff in the power spectrum mm- 
One of these theories is a kinetic-dominated (or fast-roll) 
start to slow-roll inflation, studied in [3HI2]. 

The standard slow-roll model of inflation works by 
positing a new field, (p, the inflaton, with some potential 
V ((/>). One of the most common choices of potential and 
the one we use here is U = . The field rolls down 

the inflaton potential, obeying the equation of motion 


4 > + -\- 


dVjP) 

dp 
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where H = - is the Hubble parameter and a is the scale 
factor. The second term in this equation can be viewed as 
a friction term, and so the field can approach a terminal 
velocity and the field is in slow-roll. If the field starts 
very high up on a suitably chosen potential, it will rapidly 
approach slow-roll, and then will have a long period of 
slow-roll inflation until the field exits the slow-roll regime 
as it nears the minimum of the potential. 

It is possible that the effective field theory descrip¬ 
tion of inflation breaks down near the Planck scale. In 
some models this manifests itself as a cutoff on the mini¬ 
mum value of the wavenumber k in the primordial power 
spectrum This cutoff in the primordial power spec¬ 
trum also appears in theories with a finite number of 
total e-foldings of inflation. Our work is partially moti¬ 
vated by holographic ideas that suggest an upper bound 
to the total number of e-foldings of inflation [ISHTT] , and 
in Sec. US we discuss how our model relates to those 
bounds. 

In this work, we start the motion of the inflaton field 
on the potential in a period of fast-roll before slow-roll 


inflation, and vary the initial conditions of the field. We 
solve for the predicted CMB power spectrum using a full 
numerical approach, and determine which parameters are 
the best fit to the Planck data by varying them with a 
Markov chain Monte Carlo (MCMC) sampler. 


The Planck collaboration has just published its lat¬ 
est data and results [3], and has included an analytic 
model [6] intended to be an approximation to the numer¬ 
ical model we analyze here. Our model is very similar, 
but more complete by not making as many approxima¬ 
tions, and performing the full calculations for both the 
scalar and tensor primordial power spectra. Our more 
accurate model better tests the underlying theory and 
a comparison of our model to theirs^ can be found in 
Sec. IV Unfortunately, the improvements introduced in 
our model do not yield a significant improvement in the 
fit to the Planck data. 


II. BACKGROUND 


We first consider a model of inflation with a poten¬ 
tial V = We numerically calculate the primor¬ 

dial power spectrum using the Mukhanov-Sasaki equa¬ 
tion [m [19] with a Bunch-Davies vacuum [20] at early 
times and kinetic-dominated initial conditions to start 
inflation. 

Evolution of perturbations and the primordial power 
spectrum are given by the Friedmann equations and the 
Mukhanov-Sasaki equation: 




^ The authors of Ref. [6] also perform an exact numerical calcu¬ 
lation of the scalar primordial power spectrum to fit data from 
WMAP using a grid method but this exact numerical solution is 
not used by the Planck team. 
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where K is the curvature which we set to zero to require 
a flat universe, and yk are the scalar and tensor mode 
functions respectively for a mode of wavenumber fc, and a 
prime represents the derivative with respect to conformal 
time rj, while a dot represents the derivative with respect 
to proper time t. Here we are working in units where 
h = c = 8nG = 1. 

For the quadratic inflation potential in a flat universe 
with a fast-roll start where (jP ^ mP(jP^ we choose the 
initial conditions of the background for the fast-roll start 
to inflation as in [ 6 ]. We set the scale factor a = 1 and 
the conformal time ry = 0 at the start of the numerical 
integration in the fast-roll regime: 

a = \/l + ‘2hri 
^_a^_ -h^ 

z a {l + 2hr])‘^’ 

where h is the conformal time Hubble parameter at ?; = 0. 

We assume initial conditions for the scalar and tensor 
mode functions as done in 

( 2 ) 

where Hq denotes the Hankel function of the second 
kind with index zero. These initial conditions have been 
obtained by solving Eqns. @ and (§ assuming Eq. Q 
and requiring that we must have consistency with the 
predictions of inflation without a kinetic stage in the limit 
that the kinetic stage is pushed infinitely far into the past. 

To make the calculation easier for large k, we make the 
reparameterization 

^ ^ ( 10 ) 

( 11 ) 

which results in a set of transformed Mukhanov Sasaki 
equations 


X" - 2ikX' - —X = 0 

Z 

( 12 ) 

Y” - 2ikY' - —Y = 0, 
a 

(13) 

with initial conditions 
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We choose to start the numerical integration in the ki¬ 
netic dominated regime when 00 is one hundred times 
the value of m 0 o- For convenience, in several equations 
we use Tinit = In. order to set up the numeri¬ 

cal integration, we take the limit rj ^ 0 and the initial 
conditions are 


a(0) =00 = 1 
0 ( 0 ) = 00 
0 ( 0 ) = 00 = ‘^init'^4^0 



f{x) = (1 -I- 2ix)HQ (x) — 2xWi (a:). 
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We solve the Friedmann equations to get the back¬ 
ground geometry (shown in Fig. [^, and then use the 
results of this to find z"jz^ from which we solve the 
Mukhanov-Sasaki equation where we use ordinary time, 
not conformal time, as the independent variable. Finally, 
we obtain the scalar and tensor power spectra for primor¬ 
dial perturbations: 


Vn 
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The evolution of modes of k can be seen from Fig. 
and for comparison the same plot for standard slow-roll 
inflation is shown in Fig. H When » z"/z, does 
not evolve, but z does, and is roughly proportional to a 
in both the fast-roll and slow-roll cases (though not dur¬ 
ing the transition between them), causing the downward 
slope of the large k modes on Fig. When kP ^ z" jz^ 
TZ — ^ does not evolve, and we say that the modes 
have frozen out. The integration may be stopped when 
the curvature perturbation TZ stops evolving (see Fig.[^, 
which in the standard slow-roll model is usually assumed 
to happen after k/aH 1, which is a good approxi¬ 
mation for modes larger than the cutoff. For modes k 
smaller than the cutoff, freeze-out of modes can be well 
approximated by the condition that the slow-roll param¬ 
eter obeys rj <C 1 m P For these small k modes in our 
model, we find that 77 ^ 1 when t ^ 1 /m and stop the 
integration accordingly. 

The end of inflation is assumed to occur when 1/aH 
reaches a minimum. At this time, we assume instant 


^ Equations 0 and § are correct in the r ^ oo limit. In practical 
terms, the error caused by using these equations for r = 100 is 
negligible. 

^ The general condition for freeze-out is always <C z" jz. 
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FIG. 1: Comoving Hubble length as a function of scale factor 
from the numerical code for a universe with an initial period 
of kinetic-dominated fast-roll (KD), followed by slow-roll in¬ 
flation and reheating. In this plot, modes with wavenumber k 
can be represented as horizontal lines, and a range of observ¬ 
able scales are shown. We have chosen this plot to correspond 
to the best fit parameters to the Planck data, (f>Q = 20.65, and 
m = 6 X 10“® corresponding to roughly 65 total e-foldings of 
inflation. 



FIG. 2: The power in modes of TZ (shown for the best fit val¬ 
ues of m and (jio) as a function of scale factor for selected 
modes with log^Q(fc/Mpc“^) £ {—6, —5, —4, —3, —2, —1,0} 
from lower to upper. The features near ln(a) « 2 — 3 origi¬ 
nate from the transition between the kinetic stage and slow- 
roll. We see from this plot that power is suppressed in modes 
with logjQ(fc/Mpc“^) < —3.5. We show the case of standard 
slow-roll inflation in Fig. for comparison. 


reheating into a radiation dominated universe as shown 
in Fig. [2 A longer period of reheating would impact the 
matching of the solution for the power spectrum from 
inflation onto the subsequent evolution of the universe. 
As such, it would result in a shift in the relationship 
between (j)o and the value of the cutoff in k. A longer 
period of reheating would correspond to fewer e-foldings 
of inflation, and hence a decrease in the effective value 
of Us for the part of the power spectrum with k greater 
than the cutoff.^ An example of the power spectra for 


^ Since the constraints on (po are found to be weak in Sec. |IV| a 
full analysis of the impact of uncertainties in reheating was not 
performed. 



FIG. 3: The power in modes of TZ for standard slow-roll in¬ 
flation as a function of scale factor for selected modes with 
logiQ(fc/Mpc“^) G {—6, —5, —4, —3, —2, —1,0} from lower to 
upper. The convergence of the curves to approximately the 
same value shows the approximate scale-invariance of the 
power spectrum. 


both scalars and tensors is shown in Fig. 



FIG. 4: The primordial power spectrum for scalars and ten¬ 
sors from the best ht numerical solution of the Mukhanov- 
Sasaki equation. Here = 6 x 10~® and (j)o = 20.65. 


III. NUMERICAL IMPLEMENTATION 

We solve the full Mukhanov-Sasaki equation starting 
from initial conditions, and solve for exact primordial 
power spectra for both scalars and tensors, and use this as 
numeric input for the Boltzmann code CLASS [22l [23] to 
calculate the Ce’s. We perform a full MCMC calculation, 
solving for the exact numerical solution at every choice 
of parameters. We find that the speed of the program is 
roughly halved by adding the numerical computation of 
scalar and tensor primordial power spectra. 

We split the equations for the perturbations X and Y 
into real and imaginary parts, and integrate these equa¬ 
tions together with the Friedmann equations for the evo¬ 
lution of the background simultaneously. For this pur¬ 
pose, we use a C implementation of a numerical dif- 
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ferential equation solver that automatically switches be¬ 
tween methods suitable for stiff and non-stiff equations 
(LSODA) [21] and run this solver from the command line 
in the external V{k) module |2S] of CLASS. 

A comparison of our model to the data for the CMB 
power spectrum is shown in Figs. [^and[^ 



FIG. 5: A comparison with data of different predicted CMB 
power spectra with varying cutoffs (calculated using the 
CLASS Boltzmann code). The plots are all generated with 
= 6 X 10“® and where 00 £ {20.5,20.6,20.7,20.9} (from 
the lowest curve to the highest on the left hand side). 



FIG. 6: A comparison with data of different predicted CMB 
power spectra with varying cutoffs (calculated using the 
CLASS Boltzmann code). The plots of the numerical primor¬ 
dial power spectra (solid) are all generated with m,p = 6x 10“® 
and where 0o G {20.5, 20.6, 20.7, 20.9} (from the lowest curve 
to the highest on the left hand side). The dashed curves show 
the ansatz of Ref. [6], with ln(fcc/Mpc ^) G {—7,—8,—9}. 
Though all curves are capable of suppressing the lower multi¬ 
poles, they are unable to reproduce the deficit of power seen 
in the data at ^ « 20 — 40 while preserving the good fit to the 
rest of the power spectrum. 


IV. MCMC 

We run the MCMC using the CosmoSLik sampler [26] , 
and vary over the values of m and 0Oj in addition to the 


standard ACDM parameters using the likelihoods for the 
temperature and low-£ polarization CMB power spectra 
from the Planck 2013 data. The results of the MCMC 
are shown in Fig. Because large values for 0o are 
observationally indistinguishable from standard slow-roll 
inflation and in order for the chain to converge, we im¬ 
pose an upper bound on the value of 0o shown in Fig.j^ 
We find a slight preference for a kinetic start to infla¬ 
tion with the potential as opposed to always being 

in slow-roll. However, the likelihood is only a factor of 
two larger for the best fit cutoff model, and hence the 
improvement is not significant. The Planck collabora¬ 
tion |3] showed that in a similar model |S], the expected 
improvement in log likelihood is not significant compared 
to the expected improvement from fitting an ensemble of 
cosmologies with perfect power laws as initial conditions 
and departures due solely to cosmic variance. The model 
used for the power spectrum in that case is similar to our 
exact numerical solution. Our numerical solution shows 
that there is a sharper suppression of power at small k, 
and more narrow oscillations of greater amplitude in the 
cutoff region as compared to the model of |6] used in the 
Planck paper as shown in Fig. [^ Also it appears that 
the kinetic start to inflation is unable to reproduce the 
shape of the observed dip in the CMB power spectrum 
in the neighborhood of £ « 20 — 40. Interestingly, our 
best fit value of 0o ~ 20.65 is very close to saturating the 
holographic bounds of Banks-Fischler m and de Sitter 
Equilibrium [T4UI6] shown in Fig. 

V. R + R^ INFLATION 

We also adapt the method utilized for the poten¬ 

tial to an i? -I- gravity model m- This model is the 
first inflationary model proposed, and is still an excellent 
fit to the Planck data.® It has an action 



We work in the Einstein frame where the inflationary 
potential becomes 

V(0)=A4(e-'/^'^-l)^ (28) 

This model is attractive due to its prediction of re¬ 
duced tensor power. Eor this model, we also accommo¬ 
date a variable length of reheating through varying an 
additional parameter k^ corresponding to the minimum 
value of aH which occurs during the transition from ki¬ 
netic domination to slow-roll. We incorporate the same 


® Interestingly, by putting in a kinetic transient we deviate from 
the original spirit of m- That paper offered a single “correct” 
solution for cosmology and explicitly rejected other solutions, 
such as the ones to which our transients belong, as uninteresting. 
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FIG. 7: The MCMC likelihood distributions of parameters. The mean values of parameters are displayed along with one and 
two sigma confidence intervals. The parameter 00 is not sufficiently constrained. Small values of 0o are ruled out as they would 
give too much power suppression, and large values would be observationally indistinguishable from standard slow-roll inflation 
(shown by the flat part of likelihood for 0o). One can see in the 1-D 0o plot that the likelihood for the best fit (at 0o = 20.65) 
is only a factor of two larger than for the standard slow-roll inflation. 


formalism used for the model and obtain parameter 
constraints on A, kc, and 0o as shown in Fig. 10 Because 


of the reduced tensor power, this potential is a better fit 
to the Planck data than the potential; however we 

still do not obtain a significant improvement in the fit 
due to the presence of a fast-roll start to inflation. 


VI. DISCUSSION 


Our work suggests that adding in a kinetic-dominated 
start to inflation does not significantly improve the fit 
to the CMB data from Planck. Despite this, there is a 
slight (though not significant) preference for a cutoff at 
approximately the value expected from a theory of finite 
inflation such as in [28] . 

Though our model has many similarities to the analyt- 
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FIG. 8 : The primordial power spectrum for scalars from the 
numerical solution of the Mukhanov-Sasaki equation (dotted), 
and compared to the ansatz of Contaldi et al. (lower solid) 
and the analytic model again from Contaldi et al. used in the 
Planck paper (upper solid). The dotted numerical curve has 
m 0 = 6 X 10“® and (f)o = 20.5. 


ical model of Ref. [B] used to fit the Planck data in [3] , 
we had hoped that by correcting for these differences in 
the MCMC we would improve the fit to the Planck data. 
However, it appears that the Planck data is not able to 
distinguish these differences, and thus our results do not 
signihcantly favor a cutoff in the primordial power spec¬ 
trum. 
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FIG. 9: Comoving Hubble length as a function of scale fac¬ 
tor from the numerical code evolved from the fast-roll start 
to the end of inflation joined with the analytic solution for 
instantaneous reheating into radiation domination followed 
by matter domination and cosmological constant domination 
(solid). Increasing the initial value of the scalar field causes 
slow-roll inflation to begin earlier as shown by the dotted lines 
where (j)o £ {20.0, 20.5, 21.0, 21.5}, and m = 6 x 10~®. Chang¬ 
ing the value of m in a manner consistent with the observed 
perturbation amplitude results in subtle differences that are 
too small to be seen on the scale of this plot. Modes with 
wavenumber k can be represented as horizontal lines, and the 
holographic bounds of Banks-Fischler (lower) and from de Sit¬ 
ter Equilibrium (upper) are shown as dashed lines. The solid 
line corresponds to the best ht value of 4>o ~ 20.65 from the 
MCMC, and this value almost exactly saturates the Banks- 
Fischler bound. 
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FIG. 10: The MCMC likelihood distributions of parameters for the R + n inflation model. The mean values of parameters are 
displayed along with one and two sigma confidence intervals. The parameters A and (j>o together determine both the amplitude 
of primordial perturbations and their effective spectral index. The parameter kc is not sufficiently constrained. Small values of 
kc are ruled out as they would give too much power suppression, and large values would be observationally indistinguishable 
from standard slow-roll inflation. As with rn^cf^ inflation, the presence of a cutoff can at best increase the likelihood by a small 
factor relative to standard slow-roll inflation. 























